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We find that the radial gradient of anomalous O remains constant during 1977-85 
at — 10- 15 %/A U although the intensity changes by more than a factor of 100. 

These results can be used to deduce that most of the modulation of the intensities 
of these particles is occurring beyond 27 AU. We also find evidence for a latitudi- 
nal gradient of — +3 %/degree at low energies (7.1-I0.6MeV/nuc). 

1 . Introduction. In this analysis we make use of measurements of the spectrum of O nuclei from 
the Cosmic Ray Subsystem (CRS) on the Voyager 1 (VI) and 2 (V2) spacecraft [Stone et al , 
1977] and from the Goddard-University of New Hampshire experiment on Pioneer 10 
[McDonald et al . , 1977] to study the radial and latitudinal gradients and temporal variations of 
anomalous O. Quiet-time data from 1977 to 1985 are utilized. 

2. Observations. The general features of the 1 1-year cosmic-ray modulation cycle may be exam- 
ined by using the counting rate of particles which penetrate the high-energy telescope on VI as 
shown in Fig. 1 (curve "P", mainly protons >75 MeV). A change in this rate by a factor —2 is 
observed between 1977 and the minimum in 1981-82. This change is mainly due to solar 
modulation; however, some gradient effects are present in the data as well since VI is moving 
outward from 1 to — 22 AU during this time period. Figure 1 also shows a similar plot for 40- 
106MeV/nuc C+O nuclei for 24 quiet-time intervals between 1977-85 (curve "C+O"). Here the 
change in intensity between 1977 and the minimum in 1981 is a factor —5. Finally, in Fig. 1 
the temporal variations for 5.6-17.2 MeV/nuc anomalous O are shown for the same 24 quiet- 
time intervals (curve "O"). All three components show a similar pattern of variations but for the 
anomalous O the overall intensity variation is larger than the others, approaching a factor — 100. 

In order to minimize contamination by solar and interplanetary particles and to help 
separate temporal and radial variations, six quiet-time intervals labeled A-F in Fig. 1 were select- 
ed by setting limits to the maximum low-energy helium flux, in a manner similar to that 
described in Cummings et al [1984 ]. The anomalous O spectrum is derived from the observed 
total O spectrum for a particular time interval by subtracting both a low-energy solar or inter- 
planetary component and the high-energy galactic cosmic-ray component [ Cummings et al., 
1985]. 

In Fig. 2a, b, c we show the differential intensities of anomalous O nuclei in three separate 
energy intervals for each of the six quiet-time intervals at each of the spacecraft as they moved 
outward from the sun. The solid lines are least-squares fits to the data points, except for the 
lowest energy interval for periods E and F when VI is at —24 degrees north heliographic lati- 
tude. (VI left the ecliptic plane after encounter with Saturn in 1980.) In these two cases the VI 
points, connected by the dashed lines, are significantly above the straight lines (representing a 
constant radial gradient) connecting the V2 and P10 points. We attribute the deviation to a lati- 
tude gradient at low energies having a weighted average value of 3.0 ± 1.0%/degree for the com- 
bined two periods E and F. 

In Fig. 3 we show the calculated radial gradient (slope of solid lines in Fig. 2) for each of 
the three energy intervals. In the lowest energy interval (7.1-10.6MeV/nuc) the spectrum of 
anomalous O is changing with time [Cummings et al, 1985], which may contribute to the ob- 
served variation in the radial gradient. Above 10.6 MeV/nuc no such energy spectral changes are 
found and the resulting gradient is remarkably constant. In the 10.6- 17.1 MeV/nuc interval the 
average value of the gradient is 10.7 ± 0.6%/AU, somewhat lower than the 15 ± 3%/AU found by 
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Fig. 2. Intensity of anomalous O versus heliocentric radial distance in three energy intervals from 
the VI, V2, and P10 cosmic-ray instruments. The labels A-F refer to the six quiet-time intervals 
chosen for analysis. The lines are described in the text. 
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Webber el al. (1981] for the period 1972-79. 
It is clear from Figs. 2 and 3 that the radial 
gradient for all three energy intervals has not 
changed appreciably with time or radial dis- 
tance, maintaining a value between — 10- 
15%/AU despite the fact that the intensity 
has varied by more than a factor of 100. 

3. Discussion. We believe that the constancy 
of the gradient with respect to time both be- 
fore and after the solar magnetic field rever- 
sal in 1980, and during a period when the in- 
tensity change was a factor ~ 100 is an im- 
portant clue and constraint on the solar 
modulation process. In conventional modu- 
lation theory, if the effects of particle stream- 
ing in the interplanetary medium can be 
neglected, then the radial gradient, G p and 
the radial diffusion coefficient, K r , are related 
by G r -CV/K r . The average solar wind velo- 
city, V, and K r change only slightly over a 
solar cycle (see Hedgecock [1974] and Feld- 
man et al. [1979]). Also the Compton- 
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Fig. 3. Measured radial gradients of anomalous O 
versus time for the three energy intervals of Fig. 2. 


Getting coefficient C, which is related to the spectral shape, does not change appreciably since the 
spectral shape remains similar above — 10 MeV/nuc. Thus, if the gradient can be described by 
these local parameters, it would be expected to be constant as observed. 


In this simple conventional modulation model in which K r is independent of radial dis- 
tance r, the particle intensity j at r is given by j-j b exp[CV(r-r b )/K r ], where j b is the intensity at 
the modulation boundary. Since C, V, and K r are not changing appreciably with time, and as- 
suming j b is constant, the large modulation at a given position would require a time variation in 


the boundary location r b . The magnitude of the required change in the boundary distance can 
be estimated from Fig. 2b. The dotted line is an extrapolation of the intensity in the 10.6- 
17.1 MeV/nuc energy range during the time of solar minimum to the position of P10 (27 AU) 
during the solar maximum period of interval D. The implied intensity change from period A to 
D at 27 AU is a factor of 1 36 ± 52, indicating that most of the modulation during period D is oc- 
curring beyond 27 AU. Using the average radial gradient for this energy interval we find that the 
required change in the boundary distance r b is 46 ±4 AU. A similar boundary shift has been sug- 
gested by Evenson et al. [1979] to explain electron observations. 


Such a boundary shift would also produce changes in the intensity of other particle 
species. For example, the penetrating particle rate (P) in Fig. 1 shows a variation from solar 
minimum to solar maximum of a factor —2. The median rigidity of the particles dominating 
this rate is —1.8GV, the same rigidity as 7 MeV/nuc anomalous O if O is singly charged (as ex- 
pected if they are freshly-ionized neutrals [Fisk et al., 1974]). The spectral shapes of the 
anomalous O (above lOMeV/nuc) and high-energy protons are similar (spectral index — 2) im- 
plying they have similar values of C. Since K^p^R), where p is the particle velocity and R is 
the rigidity, the expected gradient for the P rate can be scaled from the gradient of anomalous O 
by G/PJ-Gj/OXPo/Pp) — 10(.12/.89) -1.5%/AU, a value in approximate agreement with ob- 
servations by others (see McKibben et al. [1982] and Lockwood and Webber [1984]). The inten- 
sities at two different times are related by ji/j2-exp(G r (5r b )) Therefore, a boundary change of 
46 AU and a radial gradient of 1.5-3%/AU would result in j ( /j 2 -2-4 for the penetrating parti- 
cles, similar to what we observe. 


An alternative way to accomplish the same modulation without changing the boundary 
distance is to decrease K r , and thus increase the gradient, in a localized shell of turbulence in the 
outer heliosphere, as suggested by Burlaga et al. [1984]. The observed change in modulation 
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would require, for example, an increase in the anomalous O gradient from 1 1 %/AU to 44%/AU 
in a shell 1 5 AU thick. 

It is also possible that variations in the source of the anomalous O may contribute to the 
observed solar cycle variation. We note that Jokipii 11985), using a model of the anomalous 
component in which the particles are accelerated at the polar termination shock and drift to the 
solar equatorial regions, has calculated spectra of anomalous O that depend on the polarity of the 
magnetic field and which are in reasonable agreement with observations of changes in the spectra 
associated with the magnetic field reversal (see Cummings el al. [1985]). However, we find no 
evidence of a change in the gradient of the anomalous O at the time of the solar magnetic field 
reversal as would presumably be expected from such a model. 
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